Relativistic Fluctuations and Anomalous Darwin Terms in Superconductors 
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The anomalous Darwin term, one of the recently derived relativistic corrections to the conventional 
theory of superconductivity, is analysed in detail. A new derivation of this term, much simpler than 
that given originally, is presented and used to interpret the term physically, in terms of fluctuations of 
the relativistic wave packet describing the paired electrons. Several puzzling features of the original 
derivation find a simple explanation in terms of this reinterpretation. The question of observability 
of the anomalous Darwin term, and its relation to the conventional Darwin term and to earlier 
proposed superconducting Darwin terms are also discussed. 
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Recently it was demonstrated that a number of inter- 
esting and potentially observable relativistic effects exist 
in superconductors. The corresponding relativistic the- 
ory of superconductivity was proposed in Ref. de- 
veloped in detail in Refs. and ||^, and briefly sum- 
marized in Ref. Q. A recent partial review is found 
in the textbook Ref. Among the predictions of the 
theory are the existence of new types of unconventional 
superconducting order parameters (not derivable from 
Schrodinger's equation but only from the Dirac equa- 
tion), and weakly relativistic correction terms to the con- 
ventional theory of superconductivity. The present re- 
port is concerned with the physics of one of these correc- 
tion terms. 

Two essentially different types of such terms turned out 
to exist. The first class are those corrections known from 
the application of relativity to other areas of solid-state 
or atomic physic, such as the conventional kinetic-energy, 
spin-orbit, and Darwin corrections. The latter two terms 
depend crucially on the presence of an external electro- 
static potential (in the present case a crystal lattice) , and 
are zero for free particles. The second class comprises 
terms (of the same order in 1/c and formally of a sim- 
ilar structure as those in the first class), which contain 
the superconducting pair potential in place of the lattice 
potential. These terms depend crucially on the presence 
of a pair potential, and are zero for unpaired particles. 
The superconducting counterparts to the spin-orbit and 
Darwin terms were termed the anomalous spin-orbit and 
Darwin terms, respectively . 

Having derived these terms, one is now confronted with 
the task to analyze their consequences and assess their 
observability. As a first step of this program the con- 
ventional and the anomalous spin-orbit coupling have re- 
cently been investigated in detail in the framework of the 
theory of dichroism in superconductors 1^,0. It turned 
out that experimental results can be qualitatively under- 
stood on the basis of that theory . Another recent ap- 
plication of the theory is the proof by Strange |^ that the 



London equation governing the electrodynamics of simple 
superconductors (which is usually derived in a nonrela- 
tivistic way) , can be derived directly from the relativistic 
theory of superconductivity. 

The purpose of the present report is to continue with 
this program by taking a look at another term of the sec- 
ond class, namely the anomalous Darwin term. A very 
simple rederivation of that term is given, which does not 
rely on the formalism of Refs. |^ and This rederiva- 
tion provides a, previously unavailable, physical inter- 
pretation of that term, which is used to explain several 
puzzling features of the original derivation. Conceptual 
consequences and the question of observability are dis- 
cussed, as is the relation of the anomalous Darwin term 
to similar terms, proposed by other authors. 

The starting point is the standard Bogolubov-de 
Gennes equation |^,^ of the microscopic theory of su- 
perconductivity, which for spin-dependent potentials be- 
comes a 4 X 4 matrix equation. For ease of notation 
this equation is here condensed in 2 x 2 form, by intro- 
ducing the Pauli matrix Uy and the unit matrix /. For 
singlet superconductors in the absence of magnetic fields 
the equation then takes the form 



h{v)i io-yA(r, r') 
-io-yA*(r,r') ~h{v)i 

~ En 



w(r) 



(1) 



where h{r) — — ?i^V^/(2m) -I- v{r) is the single-particle 
Hamiltonian describing the normal state, and v{r) is the 
lattice potential, including a mean field ^ or Kohn-Sham 
|]lO| term arising from the Coulomb interaction. A(r, r') 
is the (generally nonlocal) pair potential of the super- 
conductor, and u„o.(r) and w„cr(r) are the particle and 
hole amplitudes, coupled by this potential. This equa- 
tion, which is usually derived in a nonrelativistic way 
1^,^ , can be identified as an approximation to lowest or- 
der in v/c (where w is a typical particle velocity and c 
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the velocity of light) to the relativistic Dirac-Bogolubov- 
de Gennes equation j^J^ , which describes pairing on the 
basis of the Dirac equation. 

To second order one finds various weakly relativistic 
correction terms to Eq. The one of interest 

in the present work is the so called anomalous Darwin 
term, which is to be added to the offdiagonal entries of 
the matrix in Eq. (|^). This term is given by 



(2) 



where x = r — r' and R = (r + r')/2 are relative and 
center-of-mass coordinates of the two electrons in the 
Cooper pair, respectively. For comparison purposes we 
also record the conventional Darwin term 



rV2z;(r), 



(3) 



which appears as a correction to the diagonal entries 
of Eq. (^). This term is usually derived by apply- 
ing the Foldy-Wouthuysen transformation or the Pauli 
approximation to the Dirac equation In the 

same way, the anomalous Darwin term (together with 
the other anomalous corrections )follows from the Dirac- 
Bogolubov-de Gennes equation 

In the conventional case a much simpler alternative 
derivation exists, which is based on the fact that a 
wave packet formed from solutions to the Dirac equa- 
tion performs a very rapid oscillatory motion with am- 
plitude Sr w Ac := Ti/(mc), where Ac is the reduced 
Compton wavelength. This so called 'Zitterbewegung' 
(trembling motion) results from the interference of the 
positive-energy and negative-energy components of the 
wave packet pJlT|-|r^. On the average, the particle thus 
feels the potential spread out over a region of order (5r, 
instead of its value at r. The resulting average change in 
energy is easily estimated to be pl]-p^ 



(5u(r)) = {v{r + (5r) — v{v)) 
dV Iv-. . \ 



\{6vfVlv{v) 



(4) 



where it was assumed that the fluctuations do not have 
a preferred direction in space. The resulting average po- 
tential {5v{t)) is, up to a numerical factor 6/8, that of 
the conventional Darwin term (r) . (The more accurate 
estimate {5y)^ ~ (3/4)Ac, occasionally found in the lit- 
erature reproduces the Darwin term exactly.) The 
above line of reasoning is very powerful, in spite of its 
simplicity. It requires, for example, only a slight general- 
ization to yield an accurate estimate for the Lamb shift 
(see, e.g., Ref. jl^, p. 80), which is a much more subtle 
relativistic effect than the Darwin term. 



We now proceed by generalizing the same procedure to 
the superconducting case. Due to the presence of the pair 
potential A(r, r') the particles now fluctuate in two po- 
tentials, each of which affects the single-particle energy. 
In a first approximation the fluctuations in the potentials 
w(r) and A(r,r') can be treated as independent. It thus 
remains to calculate the average change of energy due to 
the fluctuations in the pair potential. In the same way as 
above one finds for the fluctuations of the center-of-mass 
coordinate of the pair 

(5A(x, R)) = (A(x, R + 5R) - A(x, R)) = 
OA(x,R) l>r^,„,„ 9^A(x,R) \ 



dR^dR, 



«-(<5R)^V^«A(x,R). (5) 

Using again the estimate (5R)^ ~ for the amplitude 
of the Zitterbewegung, we find the effective potential 



((5A(x, R)) 



:V?eA(x,R), 



(6) 



which is, up to the same numerical prefactor as in the 
conventional case, the anomalous Darwin term, as ob- 
tained from the microscopic theory . This simple red- 
erivation now allows us to draw a number of conclusions 
regarding the nature of the anomalous Darwin term. 

First, note that the anomalous Darwin term is one of 
the main results of Refs. j|] and obtained there af- 
ter lengthy algebraic manipulations based on the general 
theory of Ref. ||] . The fact that one can find a much sim- 
pler rederivation of this term, not based on the general 
equations of those papers, gives additional confidence in 
that theory. 

Second, from the point of view of those references it 
remained a puzzle why no anomalous Darwin term con- 
taining derivatives with respect to the relative coordinate 
exists, in particular since there is such a term in the case 
of the anomalous spin-orbit coupling (there are, effec- 
tively, two anomalous spin-orbit terms, one containing 
V:j.A(x,R), the other VflA(x,R) [||). It was already 
verified in Ref. 1^ that the hypothetical existence of an 
anomalous Darwin term containing A(x, R) is incom- 
patible with the correct local limit, but the physical rea- 
son for the absence of this term was not apparent from 
the microscopic calculation of that work. 

The present rederivation of this term from relativis- 
tic fluctuations of the paired particles in the pair po- 
tential suggests a simple solution to this puzzle: Once 
the electrons are paired, their fluctuations are not in- 
dependent anymore. The electrons in the Cooper pair 
are locked into mutually time-conjugate single-particle 
states, so that relative fluctuations within the pair are 
strongly suppressed. We are thus dealing only with fluc- 
tuations of the center-of-mass coordinate R of the pair, 
but not of the relative coordinate x. Essentially, the en- 
tity fluctuating in the pair potential is the pair as a whole. 
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Additional insight into the nature of the anomalous 
Darwin term is gained by comparing it with the supercon- 
ducting Darwin terms derived in Refs. and In 
these papers it is pointed out that even in a nonrelativis- 
tic situation the Bogolubov-de Gennes equations do have 
negative-energy solutions (representing the condensate), 
separated by an energy gap from the positive-energy 
states (the excited quasiparticles, representing broken 
pairs). By approximately decoupling these two types 
of states via a Foldy-Wouthuysen transformation, purely 
superconducting counterparts to the Zitterbewcgung and 
the conventional Darwin term were found. These super- 
conducting Darwin terms are superficially similar to the 
one discussed here. However, although of superconduct- 
ing origin, they do not contain the pair potential, but the 
ordinary lattice potential, and are thus closer related to 
the conventional Darwin term than to the anomalous one. 
Furthermore, not being of relativistic origin, they do not 
carry the relativistic prefactor / (Sm^c^) of the Dar- 
win terms discussed here. Finally, the superconducting 
Darwin terms of Refs. and do not appear in the 
standard Bogolubov-de Gennes equation (^, but only in 
an approximate form of that equation, in which the ex- 
cited quasiparticles were approximately decoupled from 
the condensed pairs [n4 15|. This is completely analogous 



to the case of the Dirac equation, from which the conven- 
tional and the anomalous Darwin terms are derived by 
approximately decoupling electrons from positrons fl^ j. 

All of the above considerations were based on the rela- 
tivistic generalization of the BCS (singlet) order param- 
eter. However, the general theory of Refs. and Q also 
provided the form of the relativistic corrections to triplet 
superconductors, and predicted the possibility of a new 
type of order parameter, involving charge conjugation 
instead of time reversal as pairing symmetry. Recently 
these terms have been analysed in more detail . While 
for triplet superconductors one finds similar anomalous 
Darwin terms as in the singlet case, for the order pa- 
rameters involving one charge conjugation operation no 
anomalous Darwin terms are found at all, not even one 
involving the center-of-mass coordinate. 

The above rederivation of the anomalous Darwin term 
provides a simple explanation for this absence. It was 
already pointed out in Ref. Q that pairs described by an 
order parameter involving one charge conjugation oper- 
ation are intrinsically relativistic; the corresponding pair 
potential vanishes in the nonrelativistic limit. On the 
other hand, the only place in which relativity enters Eqs. 
(^ and (|^) is as the origin of the fluctuations, Sr and 
(5R. In particular, to order l/c^ or, equivalently, ((5R)^, 



the pair potential in Eqs. (|5|) and (|^) is the nonrela- 
tivistic pair potential. Consequently, for order parame- 
ters involving one charge conjugation the pair effectively 
fluctuates, to order ((5R)^, only in the crystal potential, 
but not in the pair potential. Hence there are no associ- 
ated anomalous Darwin terms for these order parameters. 
This is the result also found after considerable algebraic 
manipulations from the microscopic theory . 



We now return to the conventional singlet case and 
proceed with a tentative discussion of the possibility of 
experimental observation of the anomalous Darwin term. 
The conventional Darwin term is normally a rather small 
effect. This is not only due to the prefactor but 
also related to the fact that, unlike the spin-orbit term, 
which breaks rotational symmetry, the conventional Dar- 
win term does not break any nonrelativistic symmetries. 
In the case of the anomalous Darwin term the situation 
is similar. Although this term can break symmetries be- 
cause it contains the pair potential, which may break the 
symmetry of the crystal lattice in so called unconven- 
tional superconductors [Q, this symmetry breaking is 
of entirely nonrelativistic origin and manifests itself not 
only in the relativistic correction terms, but already in 
the nonrelativistic theory based on Eq. (|l|) . It is thus ex- 
pected that experimental observation of the anomalous 
Darwin term is much harder than that of the anomalous 
spin-orbit coupling term (which may be detected, e.g., 
by magnetooptical means |||,0|). 

An essential difference between the anomalous and the 
conventional Darwin term, which may facilitate experi- 
mental detection of the former, is the temperature de- 
pendence of the pair potential, which is absent from the 
conventional Darwin term, but intrinsic to the anoma- 
lous one. Another important difference is related to 
the derivatives of the potentials. While the expectation 
value of the conventional ionic Darwin term, originat- 
ing from the ions at positions R^, is nonzero only for 
s-states, which are nonvanishing at the position of the 
nuclei 
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= -4^^P)(r-R0, 



(7) 



the anomalous Darwin term can also influence other 
states, since the pair potential A(x, R) is not a sim- 
ple sum of Coulomb-like contributions. Any signature 
of a Darwin term for non s-states (e.g., the p-state con- 
tribution to the conduction band in simple metals, or 
the partially filled d-shells of transition metals) with the 
characteristic temperature dependence of the pair po- 
tential, would constitute an unequivocal identification 
of the anomalous Darwin term. Such an identification 
may come from high-resolution photoemission or tunnel- 
ing experiments, detecting the small temperature depen- 
dent energy shifts produced by the anomalous Darwin 
term. 

A direct comparison of the order of magnitude of the 
pair potential A (of the order of 0.1 meV in BCS su- 
perconductors and about hundred times larger in high- 
temperature superconductors) with mc^ (511 keV) is 
rather discouraging, but the fact that the anomalous Dar- 
win term affects all states involved in pairing (and not 
only the s-states) partially compensates for this difference 
in energies. Furthermore, in solid-state situations the rel- 
ativistic prefactor 1/c^ docs not always imply smallness 
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of the effects, because the large number of involved par- 
ticles can make up for this factor ||l9|] , in particular in 
superconductors, where the particles act coherently. The 
contribution of the anomalous spin-orbit coupling (which 
is of the same order in 1 /c as the anomalous Darwin term) 
to the magnetooptical response of a superconductor is, 
for example, hugely enhanced near the critical tempera- 
ture by the superconducting coherence factors ^,0]. 

Due to the second derivative of the pair potential with 
respect to the center-of-mass coordinate of the Cooper 
pair, the anomalous Darwin term is larger in situations 
in which the pair potential varies rapidly in space. Ex- 
amples are thin films, surfaces, or vortices, in particu- 
lar in systems with a short coherence length. An ideal 
system for observations would thus be a copper-oxide 
superconductor, since these simultaneously have a very 
short coherence lenght is only about 10'^ Compton 
wave lengths!), a large pair potential (evidenced by the 
large energy gap), and very heavy elements in the lattice 
(which generally favours relativistic effects). 

As long as such experiments have not been performed, 
the main results of the present paper are of a concep- 
tual nature: (i) The much simpler rederivation of one 
of the main results of Refs. ||] and ||^ gives additional 
confidence in the more complex formalism of those ref- 
erences, (ii) The interpretation in terms of rapid spa- 
tial fluctuations of the Cooper pair in the pair potential 
elucidates the physical origin of the anomalous Darwin 
term, which was not apparent from the earlier work. It 
also yields a simple explanation for the previously puz- 
zling fact that the anomalous Darwin term contains only 
derivatives with respect to the center-of-mass coordinate 
of the Cooper pair, and not also derivatives with respect 
to the relative coordinate, (iii) The same interpretation 
also explains the absence of any type of anomalous Dar- 
win terms for relativistic order parameters involving one 
charge conjugation, (iv) The temperature dependence of 
the anomalous Darwin term, together with its effect on 
other than s-states, were argued to provide a clear signa- 
ture of this term. Thin films of high-temperature super- 
conductors (or other superconductors with a short coher- 
ence length and a large pair potential) were suggested as 
ideal systems to look for observable consequences of the 
anomalous Darwin term. 
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